Conservation laws for strings in the Abelian Sandpile Model 
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The Abelian Sandpile generates complex and beautiful patterns and seems to display allometry. 
On the plane, beyond patches, patterns periodic in both dimensions, we remark the presence of 
structures periodic in one dimension, that we call strings. We classify completely their constituents 
in terms of their principal periodic vector k, that we call momentum. We derive a simple relation 
between the momentum of a string and its density of particles, E, which is reminiscent of a dispersion 
relation, E = \kf . Strings interact: they can merge and split and within these processes momentum 
is conserved, ka ~ 0. We reveal the role of the modular group SL(2, Z) behind these laws. 
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Introduction : Since the appearance of the master- 
piece by D'Arcy Thompson [l], there have been many 
attempts to understand the complexity and variety of 
shapes appearing in Nature at macroscopic scales, in 
terms of the fundamental laws which govern the dynam- 
ics at microscopic level. Because of the second law of 
thermodynamics, the necessary Self- Organization can be 
studied only in non-equilibrium statistical mechanics. 

In the context of a continuous evolution in a differen- 
tial manifold, the definition of a shape implies a bound- 
ary and thus a discontinuity. This explains why catas- 
trophe theory, the mathematical treatment of continuous 
action producing a discontinuous result, has been devel- 
oped in strict connection to the problem of Morphogen- 
esis m. More quantitative results have been obtained 
by the introduction of stochasticity, as for example in 
the diffusion- limited aggregation [3|, where self-similar 
patterns with fractal scaling dimension Q emerge which 
suggest a relation with scaling studies in non-equilibrium. 

Cellular automata, that is, dynamical systems with 
discretized time, space and internal states, were origi- 
nally introduced by Ulam and von Neumann in the 1940s, 
and then commonly used as a simplified description of 
phenomena like crystal growth, Navier-Stokes equations 
and transport processes 0. They often exhibit intrigu- 
ing patterns [y], and, in this regular discrete setting, 
shapes refer to sharply bounded regions in which peri- 
odic patterns appear. Despite very simple local evolution 
rules, very complex structures can be generated. The 
well known Conway's Game of Life can perform compu- 
tations and can even emulate an universal Turing ma- 
chine (see @ also for a historical introduction on cellular 
automata). 

In this Letter, we shall concentrate on a particularly 
simple cellular automaton, the Abelian Sandpile Model 
(ASM). Originally introduced as a statistical model of 
Self-Organized Criticality Q, because it shows scaling 
laws without any fine-tuning of an external control pa- 
rameter, it has been shown afterwards to possess remark- 



able underlying algebraic properties [3, 0| , and has been 
studied also in some deterministic approaches, exactly 
in connection with pattern formation (T^]. The ASM 
has been shown to be able to produce allometry, that 
is a growth uniform and constant in all the parts of 
a pattern as to keep the whole shape substantially un- 
changed, and thus requires some coordination and com- 
munication between different parts [lH. This is at vari- 
ance with diffusion-limited aggregation and other mod- 
els of growing objects studied in physics literature so 
far, e.g. the Eden model, KPZ deposition and invasion- 
percolation d2], which are mainly models of aggregation, 
where growth occurs by accretion on the surface of the 
object, and inner parts do not evolve significantly. 

In the sandpile, the regions of a configuration periodic 
in space, called patches, are the ingredients of pattern 



formation. In llj, a condition on the shape of patch 



interfaces has been established, and proven at a coarse- 
grained level. 

We discuss how this result is strengthened by avoid- 
ing the coarsening, and describe the emerging fine-level 
structures, including linear interfaces and rigid domain 
walls with a residual one-dimensional translational in- 
variance. These structures, that we shall call strings, 
are macroscopically extended in their periodic direction, 
while showing thickness in a full range of scales between 
the microscopic lattice spacing and the macroscopic vol- 
ume size. 

The model : While the main structural properties of the 
ASM can be discussed on arbitrary graphs [sl, for the 
subject at hand here we shall need some extra ingredients 
(among which a notion of translation), that, for the sake 
of simplicity, suggest us to concentrate on the original 
realization on the square lattice , within a rectangular 
region A E 1?. 

We write i ^ j if z and j are first neighbours. The 
configurations are vectors z = {zi}i^A € {zi is the 
number of sand-grains at vertex i). Let z = 4, the degree 
of vertices in the bulk, and say that a configuration z is 
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stable \i Zi < z for all z € A. Otherwise, it is unstable 
on a non-empty set of sites, and undergoes a relaxation 
process whose elementary steps are called topplings: if i 
is unstable, we can decrease Zi by z, and increase by 
one, for all j ^ i- The sequence of topplings needed to 
produce a stable configuration is called an avalanche. 

Avalanches always stop after a finite number of steps, 
which is to say that the diffusion is strictly dissipative. 
Indeed, the total amount of sand is preserved by top- 
plings at sites far from the boundary of A, and strictly 
decreased by topplings at boundary sites. The stable 
configuration Ti{z) obtained from the relaxation of z, is 
univocally defined, as all valid stabilizing sequences of 
topplings only differ by permutations. 

We call a stable configuration recurrent if it can be 
obtained through an avalanche involving all sites in A, 
and transient otherwise Recurrent configurations 

have a structure of Abelian group @ under the operation 
z(Bw :— TZ{z + w). We have only a partial knowledge of 
the group identity for each A (see e.g. [1, 113]; recently a 
complete characterization has been achieved for a simpli- 
fied directed lattice, the F-lattice T?]) nonetheless they 
are easily obtained on a computer and they provide a first 
example of the intriguing complex patches in which we 
are interested. The maximally- filled configuration Zmax, 
with (zniax)i = z — 1 = 3 for all i, is recurrent. More 
generally, for large A, recurrent configurations must have 
average density p{z) — |A|~^ Zi > 2 + o{l) (this bound 
is tight). So structures with density p > 2, p — 2 and 
p <2 are said respectively recurrent^ marginal and tran- 
sient. 

A patch is a region filled with a periodic pattern. The 
density p of a patch is the average of zi within a unit 
tile. Neighbouring patches may have an interface, peri- 
odic in one dimension, along a vector which is principal 
for both pat ches. Let us suppose that in a deterministic 
protocol ll| we generate a region filled with polygonal 
patches, glued together with such a kind of interfaces. 
At a vertex where I > i patches meet, label cyclically 
with a = 1, . . . ,^ these patches, call pa the correspond- 
ing densities, and 9a the angles of the interfaces between 
the patch a and a + 1 (subscripts a = £ + 1 = 1). These 
quantities are proven to satisfy the relation in Q + iQ 



[Pa+i - Pa) exp{2i9a) 



(1) 



a=l 



which has non-trivial solutions only for £ > 4 . 

A string is a one-dimensional periodic defect line, with 
periodicity integer vector k — {kx,ky) G N^, that we 
call momentum, in a background patch, periodic in both 
directions, and has fc as a periodicity vector. The back- 
ground on the two sides may possibly have a periodicity 
offset. See Fig.[Tl 

Results : In this section we report about some results 
obtained by our investigations. These facts have explicit 
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Figure 1. A string with momentum (6, 1), in a background 
pattern with periodicities V = ((2, 1), (0, 2)) . String and 
background unit cells are shown in gray. The density in the 
string tile is p = (18 ■ 3 + 8 ■ 2 + 4 ■ 1 7 ■ 0)/(6^ + 1^) = 2. 



proofs, whose details will be reported elsewhere. 

Strings in the maximally-filled background emerge in a 
simple protocol. Consider a rectangular region A, and the 
configuration Zmax- Add one grain of sand at some vertex 
j. The configuration after relaxation contains an inner 
rectangle, of strings (1, 0) and (0, 1), equidistant from the 
border of A and having j on its perimeter, and its corners 
are connected to the corners of A with strings (1,1) and 
(—1,1). There is one defect exactly at j, manifested as 
a single extra grain, w.r.t. the underlying periodic struc- 
ture. See Fig. [51 This configuration is recurrent. Now 
remove this extra grain (the configuration is now tran- 
sient), and repeat the game at some new vertex j' , say in 
the region below the inner rectangle. In the configuration 
after relaxation appear also strings (2, 1) and (—2, 1), and 
once again we have a defect at j'. This procedure can be 
iterated, and, if A is large enough, strings with higher 
and higher momenta generated. Furthermore, given that 
the unit tiles of momenta are classified (see the follow- 
ing), this protocol is completely predictable, for arbitrary 
A, through a purely geometric construction. 

In a given recurrent background, with translation vec- 
tors V = {vi,V2), one and only one string of momen- 
tum k = miVi + m2V2 = mV can be produced, if 
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Figure 2. On the left, the configuration obtained after relax- 
ation from Zmax plus an extra grain of sand exactly at the 
vertex where a defect appears. On the right, the result after 
removing the defect and the addition of one more grain. 
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gcd(mi,m2) = 1, and no strings of momentum k ex- 
ists for fc not of the form above. V and V' are equivalent 
descriptions of the background periodicity iff V — AIV, 
with M G SL(2,Z). Accordingly, m' = mM''^. And in- 
deed, sets of m G with given gcd's are the only proper 
subsets invariant under the action of SL(2,Z). The gcd 
constraint arises in the classification of the elementary 
strings, because, when d — gcd{kx,ky) > 1, the corre- 
sponding periodic ribbon is just constituted of d parallel 
strings with momentum k/d. 

The unit tile of each string, as well as of each patch, 
is symmetric under 180-degree rotations. In particular, 
momenta k and —k describe the same string. The tile of 
a string of momentum k fits within a square having k as 
one of the sides, so that each string is a row of identically 
filled squares. This is a non-empty statement: the tile 
could have required rectangular boxes of larger aspect 
ratio, and even an aspect ratio depending on momentum 
and background. 

A string of momentum k has an energy E, defined as 
the difference of sand-grains, in the framing unit box of 
side fc, w.r.t. Zmax- We have the relation E = |fe|^, or, in 
other words, the unit tile has exactly marginal density, 
p — 2, irrespectively of the density of the surrounding 
background (as seen, e.g., in Fig. [TJ. 

Two strings, respectively of momentum p and q, can 
collapse in a single one of momentum fc (see Fig. [3]). In 
this process momentum is conserved: p + q = fc. More 
precisely, the strings join together in such a way that 
the square boxes surrounding the unit cells meet at an 
extended scattering vertex, a triangle of sides of lengths 
equal to |fc|, \p\ and \q\, rotated by 90 degrees w.r.t. the 
corresponding momenta: given this geometrical construc- 
tion, momentum conservation rephrases as the oriented 
perimeter of the triangle being a closed polygonal chain. 

Local momentum conservation and the fc — fc sym- 
metry are reminiscent of equilibrium of tensions, in a 
planar network of tight material strings, from which the 
name. 

On networks, this local conservation is extended to a 
global constraint. Choose an orthogonal frame (x, t), and 
orient momenta in the direction of increasing t. Then, 
sections at fixed t are all crossed by the same total mo- 
mentum. Rigid extended domain walls between periodic 
patterns, satisfying similar local an d g lobal conserva- 
tions, appear in certain tiling models I15|, which remark- 
ably show a Yang-Baxter integrable structure, where the 
corresponding strings are usefully interpreted as world- 
lines of particles in the {x,t)-irame. Note, however, that, 
at variance with these models, in the ASM we have an 
infinite tower of excitations, for a given background, and 
infinitely many different backgrounds too. 

In the maximally-filled background, because of the {D4 
dihedral) symmetry, we can restrict without loss of gener- 
ality to study strings of momentum fc with both compo- 
nents positive. For each such fc with gcd{kx, ky) = 1, sim- 




Figure 3. A scattering involving pseudo-propagators with mo- 
menta (4,0), (2,1) and (6,1), on the background pattern of 
Fig. [1] (also symbol code is as in Fig. [!}. 

pie modularity reasonings show that there exists a unique 
ordered pair of momenta p and q, with non- negative com- 
ponents, such that p + q = k and the matrix {'^Iql) is 
in SL(2, Z). We write in this case fc ^ (p, q). For exam- 
ple, (10,3) ((7, 2), (3,1)). The endpoint of p, start- 
ing from the top-left corner of the fe framing box, is the 
(unique) lattice point which is nearest to the top-side of 
the box. This alternative definition generalizes to non- 
trivial backgrounds, and the miVi + m2V2 sublattice. 

We have found a simple algorithm to derive the string 
textures. The tile for fc ^ (p, q) is essentially composed 
of four interlaced tiles, two for p and two for q, adjacent 
to the four vertices of the fc box. This opens a problem of 
consistency for the overlapping region, of side ^ \p ^ q\, 
which is solved by the fact that, for fe {p,q), then 
p ■(— {p — q,q) if IpI > \q\, and q ^ {p,q~p) otherwise, 
again by a property of SL(2,Z) (the reader can see in 
Fig. [3] the similarity of the tile of a composed string with 
the ones of its component). By iterating this procedure, 
starting with the strings of minimum momentum, we can 
understand the textures of the full catalog of strings, and, 
in particular, this makes the protocol for producing net- 
works of strings completely predictable. 

Let us go back to the problem of £ interfaces which 
meet at a given corner, but allow now, besides interfaces 
between patches, incident strings. Following the analysis 
of [ll|, we introduce the graph-vector T = {T^}, where 
Ti is the number of topplings at i in the relaxation of the 
starting configuration, and study its characteristics in a 
region that, in the starting configuration, was uniformly 
filled with a patch. However, now we allow for toppling 
distributions which are piecewise both quadratic and lin- 
ear (the linear term was neglected in [llj, as subleading 
in the coarsening). 

For any relevant direction a, allow for a patch in- 
terface, or a string, or both. Call the difference 
for unit length (not for period), in the total number of 
grains of sand w.r.t. Zmax, due to presence of a string, 
i.e. E^""^ = i5(")/|fe*^"^ I, or the contribution from a non- 
zero impact parameter in the interface. It can be shown. 
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Figure 4. The recurrent, marginal and transient patches con- 
structed from the propagator k — (3,2) ((2, 1), (1, 1)) in 
the Zmax background, having densities p = 2 and 2 ± 1/12 
(symbol code is as in Fig. [T| . 



by reasonings as in [Uj, that the difference between the 
extrapolated toppling profile for two contiguous patches, 
at a polar coordinate (r, 9) , must be of the form 



+ r E'^°'hm(0 ~ da) 



(2) 



0(1). 



patches, respectively recurrent, marginal and transient, 
through a geometrical construction, involving p and q, 
sketched in an example in Fig. SI 

Reflection and refraction events also appear, if fc <— 
{p, q), in a single string of momentum k' — mp + q (for 
m a large integer), and in the scattering of mp + q into q 
and m parallel p strings. Indeed, as a consequence of the 
recursive construction of the string textures, the string of 
momentum k' is forced to look as a strip-shaped patch 
of TO-period width, of the marginal tile associated to p, 
crossed by a 'soft' string, that reflects twice per period 
fc', up to ultimately leaving the marginal patch, through 
a refraction, and propagates in the recurrent background. 

We shall show elsewhere that the interplay between 
strings and patches, both at the level of classification, 
and of evolution in deterministic protocols, is the key- 
ingredient to clarify allometry in pattern formation for 
the ASM, and to design new protocols in which short- 
scale defects are totally absent. The resulting structure 
is a fractal, a Sierpihski triangoloid, where the theoretical 
formula ([T]) has infinitely many distinct realizations. 



Then, by summing over a and matching separately the 
quadratic and linear terms, we conclude that, for each 9, 



ELi {Pc+i - Pa) sin^e - 9a) = 



(3) 



so that, besides the anticipated equation ([TJ for patches 
alone, that was deduced in llll, we obtain 



(4) 



which describes the string and interface-offset contribu- 
tions. 

In ([1]), the first non-trivial value for £ is 4 [ll|. In 
our generalization, 4 is the minimal value for the number 
of patches plus the number of strings, and thus includes 
new possibilities: a scattering event, with three incident 
strings in a single background, as in Fig. [31 and the case 
of two strings and two patches, producing diagrams remi- 
niscent of total reflection and refraction in optics, so that 
the specialization of ^ can be read as a Snell's law for 
ASM strings. For the case of three strings on a common 
background B, we get 



^ |fc(")|2 
a — 1 ' ' 



(5) 



which shows that momentum conservation implies a dis- 
persion relation of the form — celfc^"-*]^, and vicev- 
ersa. 

The classification of the strings preludes to a classifica- 
tion of the patches. To any string of momentum fc, uni- 
vocally decomposed as fe (p, q), we can associate three 
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